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Abstract. We consider Schrodingcr operators on [0, oo) with compactly sup- 
ported, possibly complex-valued potentials in -/^([O, oo)). It is known (at least 
in the case of a real- valued potential) that the location of eigenvalues and res- 
onances determines the potential uniquely. From the physical point of view 
one expects that large resonances are increasingly insignificant for the recon- 
struction of the potential from the data. In this paper we prove the validity 
of this statement, i.e., we show conditional stability for finite data. As a by- 
product we also obtain a uniqueness result for the inverse resonance problem 
for complex-valued potentials. 



1. Introduction 

Inverse scattering theory as well as inverse spectral theory for the Schrodingcr 
equation 

-y" + qy = ^y 

are classical subjects, its central tenets having been established some 60 years ago 
by Borg, Levinson, Gelfand and Levitan, Krein, and Marchcnkc0. During this time 
a vast body of literature on the subject has been created. A particular class of prob- 
lems, not quite so well-established, are the inverse resonance problems which are 
formulated only for a much narrower class of potentials. However, they are highly 
interesting from the point of view of applications since eigenvalues and resonances 
are directly observable in spectrometers. 

We begin by expanding somewhat on the basics of scattering theorj0. To be 
specific, the equation — y" + qy — Ay will be considered on the half-line [0, oo) 
with a Dirichlet boundary condition at zero, a case of considerable importance 
(thanks to separation of variables) for potential scattering in three dimensions with 
a spherically symmetric potential. When q is integrablc there is a unique solution of 
~ y" + QV — z 2 y which behaves asymptotically like e tzx as long as z is in the upper 
half plane. This solution is called the Jost solution and is denoted by tp{z, •)■ At 
least when q is superexponentially decaying in the mean, i.e., when J Q e rx \q(x)\dx 
is finite for all positive r, the function ip(-,x) may be extended to the complex 
plane as an entire function for any fixed x G [0, oo) (sec, for example, Lemma 3 
in [2]). If, for some C > and p > 1, the potential q satisfies |g(a;)| < C e~ x , 
the growth order of ip{-,x) is at most p/(p — 1). The function ip(-,0) is called the 
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Jost function and is of central importance in scattering theory. If zq is a zero of 
the Jost function in the upper half plane then, due to its asymptotic behavior, 
the corresponding Jost solution is an eigenfunction of the Schrodingcr operator 
associated with the eigenvalue Zq. The zeros of the Jost function in the lower half 
plane are also of physical importance. They (or rather their squares) are called 
resonances. Our interest in eigenvalues and resonances stems from the fact that 
they are fundamental objects in quantum physics with a long history, dating back to 
the early days of the theory when Weisskopf and Wigner [TH] studied the behavior 
of unstable particles. Physically, eigenvalues represent states in which the particles 
are permanently localized. Resonances, however, correspond to quasi-stationary 
(metastable) states that only exist for a finite time, proportional to the inverse of 
the imaginary part of the resonance, and have energy proportional to the real part 
of the resonance. Resonances that are close to the real axis appear as bumps in the 
scattering cross section and can be measured in the laboratory. For more details 
on resonances the reader may consult Zworski [20j . 

If the Jost function V'('jO) is °f finite growth and the location of all eigenvalues 
and all resonances is known, then Hadamard's factorization theorem implies that 
it is known up to a factor c p ^ where P is a polynomial. But the coefficients of P 
are determined also since it is known that the Jost function tends to one as z tends 
to infinity along the positive imaginary axis for any potential under consideration. 
The Jost function, in turn, determines directly the norms of the Jost solutions 
associated with eigenvalues and a quantity called the scattering phase. Marchenko's 
inverse scattering theorem states that the eigenvalues, the norming constants, and 
the scattering phase determine uniquely the potential of the Schrodinger equation, 
assuming it is real-valued and satisfies the moment condition x\q(x)\dx < oo. 
Thus, as a corollary, we have that the location of all eigenvalues and resonances 
determine a real superexponcntially decaying potential. To our best knowledge this 
observation was first publicly made by Korotyaev in 2000 (published in [8]) but 
Zworski [21j had realized (but not published) it earlier in the context of compactly 
supported even potentials on R. 

We shall call the problem of obtaining a potential from just the location of 
all eigenvalues and resonances an inverse resonance problem. Inverse resonance 
problems are eminently interesting in a practical sense since, as mentioned before, 
eigenvalues and (small) resonances are attainable in the laboratory. This is in 
contrast to the scattering phase which is not easily measured. But, admittedly, 
finding all resonances — as the theorem requires — is just as elusive a goal as 
finding the scattering phase. Therefore a fundamental question arises: What may 
be said about a potential when only the location of the associated eigenvalues and 
small resonances is known. 

The following argument shows that large resonances carry little physical infor- 
mation. If the potential q is compactly supported in [0, 6], absolutely continuous on 
its support, and if q(b) ^= then it is well known that there are only finitely many 
eigenvalues and that the (roots of the) resonances x n + iy n are asymptotically close 
to the curve given by 



nutest change of the potential near the right endpoint of its support and suggests 




This shows that the asymptotic distribution of resonances changes upon the mi- 
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that one might be able to say a good deal about the potential from knowing the 
location of eigenvalues and resonances of modest size and, in particular, without 
knowledge of the asymptotic distribution of the large ones. 

Indeed, despite the fact that the finite data inverse resonance problem (or the 
finite data eigenvalue problem, on compact intervals) is ill-posed, having no unique 
solution, many numerical methods have been developed for its solution. Without a 
claim to completeness we mention here Andrew [1], Brown, Samko, Knowles, and 
Marietta [3], Hald [BJ, Paine [TB], Rohrl [T7], Rundell and Sacks [IS]. The question 
now is how close these solutions are to each other and, more importantly, how close 
they are to the actual potential. The usual way to answer this question in the 
numerical analysis literature is to apply the recovery algorithm to a situation in 
which finite spectral data were generated from some known potential. The quality 
of the recovery procedure is then assessed according to how closely the recovered 
potential approximates the original one in some norm (sometimes the "eyeball" 
norm). Mathematically this involves are large leap of faith (practically such a leap 
of faith may, of course, be necessary). 

Surprisingly it seems that this stability problem has not received much attention 
(one exception is Hitrik [7]). Even for the much simpler inverse eigenvalue problem 
on a compact interval, this question was only answered as recently as [12j . In [13] 
it was addressed for a discrete Schrodinger equation. Some stability results for the 
case of a real- valued potential are given by Korotyaev [9] but these do not address 
the case of finite resonance data, and indeed require quite delicate knowledge of 
the large-resonance asymptotics which will certainly not be available if only finitely 
many resonances are known. 

In this paper we allow the potentials to be complex-valued but assume that they 
lie in a ball of fixed radius in L l (0, 1) and are compactly supported. We suppose 
that two potentials q and q are both known to have compact support in some fixed 
interval - without loss of generality, the interval [0, 1] - and that for some R > and 
for some e > 0, their resonances and eigenvalues are e-close inside the disc of radius 
R centered at zero. Outside the disc of radius R, the resonances and eigenvalues of 
q and q need not be close at all; no assumption on the resonances outside the disc 
is made. If we also assume that q — q is in a ball of fixed radius in L p (0, 1), p > 1, 
we obtain an estimate on 



in terms of e and R which tends to zero as 1/R and e tend to zero. Even without 
the assumption that q — q € L p (0, 1) we can show that (JT]) tends to zero as 1/R 
and e tend to zero but we have no control over the rate of convergence. Note that 
in the case where all eigenvalues and resonances for q coincide with those for q we 
obtain a uniqueness result, i.e., q and q are then equal almost everywhere. To the 
best of our knowledge this uniqueness result is new for complex- valued potentials. 

Our bound is obtained using transformation operators, Hadamard factorization, 
estimates developed from Jensen's formula and some elementary facts about Fourier 
transforms. The bulk of the work, however, lies in the estimation of the transfor- 
mation operators from the resonance data. This is based on solving a non-standard 
boundary value problem for the hyperbolic PDE satisfied by the transformation 
kernels, with data given on a characteristic line. Our main result is Theorem 16.11 
together with its corollaries 16.21 and 16.31 
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2. Transformation operators 

Definition 1. Let Q > 0. By B(Q) we denote the set of all (complex valued) 
functions q G L 1 [0, oo) which have compact support in [0,1] and are such that 

Jo M < Q- 

Throughout this work we shall consider Schrodingcr operators on the half-line 
[0, oo) in which the potentials lie in B{Q). Given any q G B{Q) and any z G C we 
consider the Schrodinger equation 

— u" + q(x)u = z 2 u, x > 0. 

Since q is compactly supported in [0, 1], for each z G C this equation has a unique 
solution u satisfying the condition 

u(z,x) = exp(izx), x > 1. 

This solution is called the Jost solution and, for each fixed x > 0, it is an entire 
function of z of growth order one. We shall denote it by ip(z,x). Note that if 
Im(z) > then ip(z, •) G £ 2 [0, oo), and so if ip(z, •) happens to satisfy the Dirichlct 
boundary condition ijj(z,0) = then z 2 will be an eigenvalue of the Dirichlet 
Schrodinger operator H (q) defined by 

Ho(q)u = —u" + qu 

on the domain 

V(H Q (q)) = {u G L 2 [0,oo) | — u" + qu G £ 2 [0, oo), u(0) = 0}. 

If, on the other hand, ip(z,0) — and Im(z) < 0, then ip(z, ■) i 2 [0,oo) and so 
ip{z, ■) cannot be an cigenfunction of Ho(q). In this case z 2 is called a resonance of 
Ho(q). Thus we have the following dichotomy of the zeros of z i— > ^(z,0): 

• If Im(z) > and -0(z, 0) = then z 2 is an eigenvalue of Ho(q) with eigen- 
function ip{z, ■). 

• If Im(z) < and tp(z,0) = then z 2 is a resonance of H (q) with wave 
function - i/ , ( z : ')• 

Where no confusion will result, we shall abbreviate tjj(z,0) simply to ip{z) and call 
this the Jost function. 

A remarkable fact about Schrodingcr equations - see, e.g., Levitan [10] - is the 
existence of an integral operator K, not depending on z, which maps Jost solutions 
for one potential to Jost solutions for a different potential. Given two potentials q\ 
and q2, let ipj(z, x) be the Jost solution of the equation — u" + qju = z 2 u, j = 1, 2. 
Then there exists a kernel K such that 

POO 

ip2(z, x) = ipi(z, x) + / K(x,t)ipi(z,t)dt. 

J x 

Throughout the later sections of this article we shall require estimates of kernels 
such as K. Since, in our situation, 

,(i-x)/2 

/ \q 2 (a-P)- qi (a + p)\d(3<2Q 

Jo 

one obtains (see, e.g., Theorem 3 and Lemma 1 of [4]) that 

oo 

K(x, t) = £ K n (x, t) =: K (x, t) + H(x, t) (2) 

n=0 
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where 

1 f 1 

K o(x,t) = - (q 2 (s) - qi(s))ds 

1 J(t+x)/2 

and, for neN, 



"1 f(t~x)/2 

K n (x,t) = I / {q 2 (a-/3)-qi(a + l3))K n -i(a-l3,a + (3)d/3da 



(t+x)/2 JO 

so that 

|^(^t)|< 1(2Qr 



t + of 


n „1 






2 


+ J(t+x)/2 



\Q2{s) - qi{s)\ds 



2 n! 
and 

1 Z" 1 

\K(x,t)\<- / | ft («)- ? i(s)|dflexp(2Q[l-(t + a?)/2] + ). 

^ J(t+x)/2 

hi particular, 

\K(x,t)\<Qc^ (3) 

and K n (x,t) = K(x,t) = if t + x > 2. Notice also that K n (x,x) = for n > 1 
and hence 



if (a:, a:) = K (x,x) = - I {q 2 (s) - qi(s))ds, (4) 



2 , 

which allows the difference of the potentials to be recovered from the transformation 
kernel K. 

It is also shown in [4] that the function t i— > H t (Q,t) is absolutely continuous, 
and 

\H t (0,t)\ <CQe 2 ®, (5) 

for some constant C independent of q\ and q 2 . 

As explained in the introduction, we consider in this paper the problem of esti- 
mating the difference between two potentials q and q whose resonances are close to 
each other, if they are not far from the origin. In order to do this we adopt some no- 
tation for specific transformation operator kernels corresponding to different choices 
of qi and q 2 above: 

• K q for the transformation from a potential to a potential q. 

• L q for the transformation from a potential q to the potential 0. 

• Kq for the transformation from a potential to a potential q. 

• B for the transformation from a potential q to a potential q. 

We denote the Jost solution for the potential q by ij; and for q by ijj. Thus we have 

p2-x 

ip{z, x) — cxp(izx) + / K q (x,t) exp(izt)dt (6) 

J x 

and 

r-2-X 

ip(z, x) = exp(izx) + / Kq(x,t) exp(izt)dt. (7) 

J X 

Correspondingly, the kernel L q maps Jost solutions for potential q back to solutions 
of the free problem with potential 0: 

p2-x 

exp(izx) — tp(z, x) + / L q (x,t)ip(z,t)dt. (8) 
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Into the right hand side of j7|) we insert the expression for exp(izx) from (|8|) to 
obtain 

p2~x p2-x 

ip(z,x) = ip(z,x) + / L q (x,t)ip(z,t)dt + / Kq{x,t)ij}{z 1 t)dt 

J x J X 

2-x r-2-s 

ds / Kq(x, s)L q (s,t)tp(z,t)dt. 

J x J s 

This means that 

/•2-x 

ip(z, x) = tjj{z,x) + / B(x,t)tp(z,t)dt, 

J x 

in which 

B(x,t) = Kq(x,t) + L q (x, t) + / Kq(x,s)L q (s,t)ds. 

J x 

This expression is standard and may be found in [10] . In the special case in which 
q = q we know that B must be zero and K q must be ifg; this yields 

= K q (x,t) +Lq(x,t) + [ Kq{x,s)L q {s,t)ds. 



In particular, this gives 



B(0,t) = Kq(0,t)-K q (0,t)+ / (^(0,s)-K,(0,s))L g (a,t)ds. (9) 



We know from (|3|) that the sup-norm of transformation kernel L q is bounded by a 
constant which depends only on \\q\\i. Thus we obtain a bound on B(0, i) from one 
on (Kq — K q )(0, •). This estimate on B(0, t), in turn, will eventually yield a bound 
for B(x,t) by an iterative procedure; the difference q — q is then found from ((¥]), 
which yields 

2i?(.T,x) = / (q-q). (10) 



In order to find a bound on (Kq — K q )(0, •) we observe that, setting x = in 
eqns. ©, ® and inverting the Fourier transform, 



(Kq-K q )(0,t) = ±- I $(z)-i>(z))exp(-izt)dz. (11) 

Z7T 



Our first task, therefore, is to estimate ip( z ) ~~ f° r rca l z - We shall start by 
doing this in the case when ip and ^ have exactly the same zeros in some large 
disc; the resulting bound on \(Kq — K q )(Q, -)| is in Theorem 14.11 The case where 
the zeros inside the disc are perturbed is handled in Section [6] 

3. Estimates on the difference of Jost functions having the same 
zeros in a disc of radius r 

In this section we will derive a pointwise bound for the difference of two Jost 
functions in the interval [ — i? 1 / 3 , i? 1 / 3 ] under the assumption that eigenvalues and 
resonances in the disc of radius R coincide. More precisely, we will prove the 
following theorem. 
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Theorem 3.1. For any positive number Q there are numbers C > and Ro > e so 
that the following statement is true for any R > Rq. If q and q are two potentials 
in B(Q) and if the zeros of the associated Jost functions ip and ip coincide in the 
disc \z\ < R then 

\ip(z) -i){z)\ < CR- 1/3 
for all z satisfying — i? 1 / 3 < z < R 1 / 3 . 

This theorem will be proved at the end of the section after several lemmas have 
been established. The key is Hadamard's factorization theorem which says that 

oo 

iP(z)=z n ° c g M l[E(z/z n ) 

where uq is a nonnegative integer, g is a polynomial of degree at most one, E(w) = 
(1 — w) e w , and the z n are nonzero complex numbers. We introduce the abbreviation 

oo 

U(R,z)= [] E(z/z n ). 

\z n \>R 

We begin by establishing a preliminary estimate for the Jost function. 

Lemma 3.2. For every positive number Q there is a positive constant k such 
that the Jost function associated with any potential q G B(Q) has the following 
properties. 

(1) \4>(z)\ < * for all zeR. 

(2) \i>(z)\ < kc 2 ^ for all zeC. 

(3) If p > then \ip(z) — 1| < n/p for all z in the disc {z : \z — 3ip| < p}. 

Proof. The representation ^{z) = 1 + J Q 2 K q (0,t) e lzt dt gives immediately that 
\ip(z) - 1| < 2||A" 9 (0, •)||o cxp(2|Im(z)|) for all z e C. This proves the first two 
statements. If Im(z) > we may estimate \tp(z)-l\ by 11^(0, -)||oo e~ tIm ^ dt < 
\\Kq(0, -)||oo/ Im(2). Since the disc {z : \z— 3ip\ < p] is contained in a sector where 
Im(z) > |z|/2 we obtain the third statement. □ 

We now assume that p > 2n and introduce the function N(r) which counts the 
number of zeros of ip contained in the disc {z : \z — 3ip\ < r}. Note first that 
iV(0) = since |"0(3ip)| > 1/2. Since ip has growth order one the counting function 
can only grow linearly. In fact, Jensen's formula 

r er N(t) i f 27r 

/ -^r ± dt=— log|V(3zp + ere It )|A-log|^(3zp)|, 
Jo t 2tt J q 

the inequality N(r) < f° r t^ 1 N(t)dt, and part (2) of Lemma 13. 21 give 

N{r) < log(2K) + 6p+ 2er. (12) 

The elementary factor E{w) = (1 — w) e w satisfies | \ogE(w)\ < 2\w\ 2 as long as 
\w\ < 1/2. Therefore, thinking of w as z/z n , we are interested in an estimate on 



\z n \>B 



Zr,\- 2 . 
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It will be convenient to assume that R > 9p > 18k. Since \z n \ > R we get 
\z n \ > 3\z n — 3ip\/4 so that 

|,„-3*p|>2fl/3 J2R / 3 1 J2R / 3 1 

Using now inequality (|12[) and log(2«;) < p/2 gives S < 36/R. With the aid of the 
inequality | e" — 1| < \u\ e' u we arrive at the following lemma. 

Lemma 3.3. Let z n , n € N, denote the nonzero zeros of the Jost function ip and 
assume that R is a positive number which exceeds 18k, where k is the quantity given 
in Lemma \3/iH Then 

|n( J R,z)-l|<-|^exp(72|z| 2 /i?). 

provided that \z\ < R/2. 

Now wc return to the case of two potentials q and q in B(Q). Let k be the 
number associated to Q according to Lemma 13.21 Since we assume the zeros of ip 
and tp within the disc of radius R to coincide we get 

i>(z) Il(R,z) 
and we need to estimate exp(g — g). 

Lemma 3.4. There are positive constants Ro and c depending only on k such that 

| eSW-SW _x| < cR-i/3 

provided that R> Rq and \z\ < i? 1 / 3 . 
Proof. Suppose \z — 3ip\ < p. Since 

effW-SW _i = ( jM + fi(R,z)-U(R,z) 

n(R,z) \t/j( z ) J U(R,z) 

we get, when p 2 /R is sufficiently small, 

i e s(*)-s(*) _!| < ^! + Ef^_ 
p R 

from (3) of Lemma 13.21 and Lemma 13.31 when A and B denote suitable numerical 
constants. The two contributions to the error are in balance when we choose p on 
the order of i? -1 / 3 . Specifically, there are positive constants i?o and c depending 
only on k such that 

| e <?(*)-§(z) _!| < _£._£- 1/3 

if p 3 = i? > R and |z - 3ip| < p. 

Suppose f(z) = cxp(aiz + ao) — 1. It is easy to show that 

\f(z)\<^exp(5e/(l-e)) 

in the disc \z\ < p if |/(z)| < e < 1 in the disc |z — 3ip| < p. Applying this to the 
case at hand gives the stated estimate after possibly increasing Rq to ensure that 
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Proof of Theorem[3Ji Suppose -R 1 / 3 < z < R 1 ' 3 . By part (1) of Lemma O 



\ljj{z)-${ Z )\<K 



1 



TL(R,z) 



U(R, 



TL(R,z) 



U(R,z) 



- 1 



Using the estimates obtained in Lemma l3.3l and Lemma l3. 41 establishes the theorem 
for a C depending only on c, Rq, k, and numerical constants and hence only on 
Q. □ 



4. Large z asymptotics of the Jost functions: further results and 
consequences for transformation kernels 

In this section we assume that q and q lie in B(Q) and that q — q lies in L p [0, 1] 
for some p > 1. If p > 2 then q — q is still in L 2 ([0, 1]) and therefore we assume 
henceforth that p e (1, 2]. 

We start with dTTJ) from Section [2 



1 

2^ 



— ip{z)) exp(—izi)dz. 



K q (0,t) - K q (0,t) 
In particular, therefore, 

R 1/e 

Kq(0,i)-K q (0,t) = — ($-ip)(z)exv{-izt)dz+-^- I {ip-ip){z) exp(-izt)dz 



1 



2lT J — Rl/6 



2tt 



(13) 



The first term on the right hand side of (fT3|) will be handled using Theorem 13.1 
which yields 



— / (tp — 4>){z) exp(— izt)dz 



2tt 



The second term, 



E R {t) 



R 1 / 6 
1 

2^ 



< 



C 
R^' 



(-0 — ip)(z) exjp(—izt)dz, 



(14) 



(15) 



will be handled using asymptotics which refine the results in Lemma 13.21 and which 
we develop using the transformation equation ([6]) (with x = 0) and integration by 
parts. Following the notation in ([2]), we obtain, after an integration by parts, 



1>{z) 



iK q (0,0) i 
z 4z 



g(t) cxp(izt)dt, 



where g(t) — q(t/2) — 4iJ t (0, t). We can also write 

»if 9 (0,0) i 
V-(z) = 1 + —g(z) 

where g is the Fourier transform of g. This immediately yields, in an obvious 
notation, 



m - m = - z (K q - K q )(o, o) - ±(fl)(z). 



(16) 



Since H t (0,t) is continuous and bounded by the bound given in ((5|) g — g is in 

LP([0, 2]) where p e (1, 2] so that (if^#) is in Lp/^ 1 ) (R) by the Hausdorff- Young- 
inequality. 
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Now we substitute (fT6]) into the right hand side of IT5|) to obtain 

E R (t) = ^-(K q - K q )(0,0) f -exp(-izt)dz 
2m J\ z \>Ry 6 z 



z\>R 

1 



-(9 ~ g)(z)exp(-izt)dz. (17) 

8?r J\z\>RV e z 

The first integral in (fTTj) can be rewritten by a change of variable £ = to as 
Jj^i >tR i/ 8 | exp(— i£,)d^, and hence an integration by parts yields 

„ cos^i? 1 / 6 ) /• , 
-2* V + z / exp-«£ -i, 

which is (^((ii? 1 / 6 ) -1 ) when ii? 1 / 6 is large. When ii? 1 / 6 is small the integral can 
be estimated by taking the Cauchy principal value. In either case, there exists a 
numerical constant C± such that 



if 1 

/ — exp(— iztjdz 

J\z\>R 1 / 6 z 



2n 



The remaining integral in (jTTJ) is estimated by the inequalities of Holder and 
Hausdorff- Young so that 



— (g — g)(z) cxp(—izt)dz 

z l >R l/6 Z 



< Czip-iyVPRV-pyWtl + \\q-q\\ p ) (19) 



where the constant C2 depends only on Q by ([5]). Note that this becomes unbounded 
as p tends to one. 

Combining the estimates (TT4"1) . (TT5)) and (TIT)]) yields: 

Theorem 4.1. Suppose that q and q lie in B(Q) and that q — q is in L p [0, 1], 
p £ (1)2]; suppose that Rq is as in Theorem \3.1\ and that R > Rq. Then there 
exists a constant C (possible larger than the one used in Theorem \3.1\) depending 
only on Q such that 



(20) 



\K q (0,t) ~ K q (0,t)\ < C( P - 1)"Vp(i + \\ q - q\\ p ) min (l, -L 

where v = (p — l)/(6p). 

The relationship ((9|) gives immediately the following corollary. 
Corollary 4.2. Under the assumptions of Theorem \4.1\ we have the estimate 

log_R N 
tR v 

where v = (p — 1) / (6p) and C is a constant depending only on Q. 

Remark 1. If we do not require a rate-of-convergence estimate then we can assume 
simply that q and q lie in L^O, 1]. In order to see this, consider the term 

Fn(t) := / cxp(— izt)dz. 



15(0, t)\ < C'ip - 1)"Vp(1 + || ? - q\\ p ) m in ( 



1. 



Using the definition of the Fourier transform and Fubini's theorem, we get 

f 1 , J f n sin(z(s -i)) , \ , 
F R {t) = 2i lim / q(s) / — ^ -dz ds. 

n ^°°Jo KJr 1 / 6 z J 
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The inner integral is bounded as function of n, R, s, and t. Hence, applying the 
dominated convergence theorem twice shows that Fn(t) tends to zero as R tends 
to infinity. 

5. Estimating the difference of two potentials from the difference 

of two jost functions 

Eqn. (|10|) will yield bounds on J^.(q — q) from a bound on B(x,x) and hence, 
in particular, from a bound on B(x,t). Corollary [42] gives a bound on B(0,t). In 
order to determine a bound on B(x,t) from the bound on B(0, t) we first observe 
that, following the material in Section [21 

B(x, t) = for x + t > 2. 

In particular, 

B(x,2-x) = 0, 0<a-<l. (21) 
We shall show that this condition, combined with the knowledge of B(0,t) for 
< t < 2, determines B completely. In order to do this we derive a second integral 
equation for B. 

Following Levitan [10] we observe that the function 

^ i*x pt—u+x 

I du , 

t-\-U — X 



z (x,t) = — J du I F(u,v)dv 



satisfies the inhomogeneous wave equation 

zo,xx - z .tt = F(x,t), 

with boundary condition zo(0,t) = 0. Putting g(x) = zq(x, 2 — x) we observe that 

, , / t + x\ ( t - x 
vo{x,t) = g II — I - g I 1 — 

satisfies the homogeneous wave equation together with the boundary conditions 
vo(0,i) = 0, v${x, 2 — x) = —g{x) = —zq(x,2 — x). Hence the function u${x,t) = 
Zo(x,t) +vo(x,t) will satisfy the inhomogeneous wave equation 

U ,xx - u 0t tt = F(x,t) 

with homogeneous boundary conditions 

ii (0,i) = 0, u (x,2-x) = 0. 

Now uq(x, t) is expressed as a sum of three integrals: 

u (x,i) = h(x,t) + h(x,t) - I 3 {x,t), 

where 



1 



I±(x, t) = Zo(x,t) = — j du j F(u,v)dv, 



t — u+x 



t + X\ 1 



t+u — X 
l-(t+x)/2 r 2-u 



h(x,t)=g[l — j = 2i du I F(u,v)dv, 



t-x\ 1 



Jt+x+u 
l-(t-x)/2 i-2-u 

h{x,t)=g\ l-^r 1 ) =^ I du F(u,v)dv. 

\ * J £ Jo Jt-x+u 

These integrals all have the same integrand and elementary calculations show that 
the regions over which integrals 1\ and I2 take place are disjoint and are sub- 
domains of the region over which integral 23 takes place. Subtracting I\ and h from 
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I3 leaves an integral over the rectangle with corners (x,t), (l — (t—x)/2, l+(t— x)/2), 
(0, t + x) and (1 — (t + x)/2, 1 + (t + x)/2). We make the change of variables 

v = a + (3, u = a — f3, dudv = 2dadj3, 

and obtain 

fl r{t+x)/2 



uq{x, t) = — I da I F(a — f3, a + /3)d/3. 

J{t+x)/2 J(t-x)/2 

If we now ask for the solution w of the wave equation 

w xx - w tt = F(x,t), w(0,t) = B(0,t), w(x,2 — x) = 0, 

then, since B(0, t) = for t > 2, w will be given by 

w(x,t)=u o (x,t)+B(0,x + t). (22) 

Now the transformation kernel B is required to satisfy the wave equation 

(B x - B t ) x + {B x - Bt) t = (q(x) - q{t))B 

subject to the condition (|2T|) . and with B(0, t) known. In view of the expression 
(122p this means we should choose B as the solution of the integral equation 

B(x, t) = B(0, x + t)+ da (q(a + /3)- q{a - fi))B{a ~(3,a + f3)d(3. 

J(t+ X )/2 J(t- X )/2 

Iteration (cf. Lemma 15.11 below) shows that this solution is given by 

00 

B(x,t) = J2 B n(x,t) 
n=0 

where 

B (x,t) = B(0,x + t) 

and 

/■l r(t+x)/2 

B n+1 (x,t)= da (q(a + f3)-q(a-p))B n (a-f3,a + p)dp. (23) 

J(t+ X )/2 J{t- X )/2 

Lemma 5.1. Suppose that q and q are in B(Q) and that there exist constants Co, 
C\ and i?2 > e such that for all t G (0, 2], 

\B(0,t)\ <d+C Q min (l, -^j . (24) 

Then 

whenever ntfj and < x + t < 2. 

Proof. The proof is by induction. We first check that the estimate holds for n = 1. 
We have 

r{t+x)/2 

\q{a + /3)- q{a - f3)\\B(0, 2a)\df3 

l{t+ X )/2 J(t- X )/2 



\B!(x,t)\ < da 

J(t+ X )/2 J(t- 

< 2Q [ \B(0,2a)\da 
Jo 



Ci + Co log|^) 
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which establishes the result for n = 1. 

Next we substitute the estimate (|25[) into the right hand side of (|23|) and try to 
recover the appropriate estimate for B n +\. Since (|25p holds we have 

| B „ +lfet) |<f Cl + c„!^)W)^/' d-^.*, 

V H 2 J (n - 1)! J( t +x)/2 

which yields the required estimate. □ 



Lemma 5.2. Under the hypotheses of Lemma \5.1\ the estimate 



\B{x,t)\ < Ci + Cb 



log(R 2 ) 



(1 



(26) 



(x + t)R 2 

holds for all (x, t) in the triangle bounded by the lines x = 0, x = t and x + t = 2. 

Proof. This is an immediate consequence of the fact that B(x,t) = B(0,x + t) + 
J2n=i B n(x, t) together with the bounds (|24I25|) . □ 

Theorem 5.3. Let Qi and Q p be positive numbers and p S (1,2]. Then there 
is a positive number C , depending only on Qi and Q p , and a positive number Rq, 
depending on Qi, Q p , andp, so that the following statement is true for any R J> Rq . 
If q and q are two potentials in B{Q\) such that \\q — q\\ p < Q p and for which the 
zeros of the corresponding Jost functions are identical in a disc of radius R, then 



sup 

x6[0,l] 



(«-«) 



< C(g 1 ,Q p )(logii)( 2p - 2) /( 2p - 1 ^-( p - 1 ) 2 / (6p(2p - 1 ». (27) 



Proof. Let 7 = (p — l)/p = d>v and < 77 < 1. Lemma [5.21 with C\ = and Co 
determined by Corollary |4.2i implies the existence of a constant C 2 , depending only 
on Qi, such that 

2\B M \<C 2 ( P -l)-^a + Q P )^^:=^ 

provided that R is at least as large as Ro given by Theorem 13.11 By possibly 
enlarging Ro we have M < jQ p when R > Rq. (If M is not much smaller than 
7<5 P our final estimate will be worse or not much better than the trivial estimate 

£(?-<?) I <Q P )- 

Thus, if < n < 1 then, for all x £ [0, 1], 



(q-q) 



< 



\q-q\ + 



(q-q) 



< Q P v 7 + 



(28) 



according to equation (flT))) . Substituting 77 = (M/(7Q p )) 1 ^ 7+1 ^ < 1, the point 
where the best estimate occurs, into (|2"8")l gives the desired result. Note that 
C(Qi,Q p ) may be chosen, independently of p. □ 

Remark 2. The proof of Theorem E3] reveals that C(Qi, Q p ) = 0{Q p p /{2p ~ 1] ) as 
Q p tends to zero. 
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6. Further errors from perturbation of the resonances 

In the previous sections we considered a change of potential from q to q which 
preserved all zeros of the Jost function in a disc of radius R. In the current section 
we allow the potentials q and q to have different zeros z n and z n inside the disc of 
radius R, satisfying a uniform bound: 

\z n -z„\<e, n = 1,2,... N :=N(R). 

Since zeros of the Jost function in the upper half plane lie in a half disc centered 
at zero whose radius is bounded by a constant times Q we may assume without 
loss of generality that none of the zeros z n or z n lie on the real axis. If this is not 
true then the contours taken below for the inversions of the Fourier transforms may 
be deformed around the zeros. We therefore assume that at every point z on our 
inversion contour, and for all n, 

\z - z„\ > 1, \z- z„\ > 1. 

Define a function 

n— 1 v ' 

Then the Jost functions ip(z) and ip(z) satisfy 

«*> - "'^^fki 

in the notation of Section [3l We still have the equation 

K q (0,t) - K 9 (0,t) = f (i'(z) - i>(z))exv{-izt)dz =: I R (t) + E R (t), 
in which 

E R (t) -- 

We also still have the estimate on E R from fTS]) , (|19p giving the existence of a 
constant C\ depending only on Q such that 

\E R {t)\ ^da + llg-gymm^-L 

in which v = (p — l)/(6p). We therefore turn to estimating lR,{t), which we first 
write as 

I R (t) = - j^Wz) \l - + ^z)(W(z) - l))eM-izt)dz. 

Estimating the first term using Lemma 13.21 (part (1)) to bound ^{z), and Lemmas 
to bound 1 — »„ f -i^S p :Z j . we obtain 

\I R (t)\ < C 2 R- 1 ' 6 + 1- [ $(z)\\W{z) - l\dz 







— ip(z)) exp(—izt)dz 




Iz^R 1 ' 8 






rR 1 ' 6 






/ Mz) - 


- ip( z )) exp(—izt)dz, 


2tt. 


/-HV6 



2tt 
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for some constant C2 depending only on Q. We now apply the inequality \W — 1| < 
I log W\ exp(| log W\) together with the bound |'0( z )l < K from Lemma [3.21 fpart 
(1)) to obtain 



\I R (t)\ <C 2 R- 1/6 + k 



,1/6 



\logW(z)\faq?{\logW(z)\)dz. 



Since | log(l + x)| < — log(l — \x\) < 2\x\ whenever |x| < 3/4 we find 

N(R) 

<*E 



N(R) 

\ogW(z)\< Y, 



log 1 



-n 



(29) 



(30) 



provided the summands on the right are bounded by 3/4. Denoting the smallest 
integer which is at least as large as x by [x] we distinguish now the cases n < 
^ei? 1 / 6 ] and n > ^ei? 1 / 6 ] . (H ere we assume that N(R) > [8ei? 1 /6] ) 

since the 

case N(R) < [Sei? 1 / 6 ] will give the same bounds a fortiori.) In the former case 
we use the assumption | z — z n | > 1 to estimate a summand in (|30[) by e which we 
assume to be less than 3/4. In the latter case we observe that thanks to eqn. (fl"2"l) 
with p = 2k, we have 

\z - z n \ > \z n \ - \z\ >^-n- R 1 / 6 - 12k - Iog(2«) > -Ui 
ze 4e 

if we require i? 1 / 6 > 12k + log(2«;). Hence, if n > [Sei? 1 / 6 ] a summand in (|3"0|) is 
estimated by 4ee/n. Using (fT2)) again we find 

\logW{z)\ < 16cei? 1/6 + 8celog(i?)+e< 17cei? 1/6 . 

Using this estimate for the argument of the exponential function in (|29[) we 
obtain 

\lR{t)\ < C 2 ir 1/6 + Kexp(17eei? 1/6 ) / \\ogW{z)\dz. (31) 
To estimate J \ \ogW\dz we stick to the bound in (|30[) . i.e., 

\\ogW{z)\dz<2e{S l + S 2 ) 1 



where 



R 1 / 6 



Si 



ISeR 1 ' 

E 

n=l 
N(R) 

E 

= r8eHV6]+l 



dz 



-Ri/e \z — Z n \ 



-Rl/6 \Z ~ Z n \ 



In the sum Si we set £ = Re(z — z n ) and observe that \z — z n \ > yl + £ 2 with 
|£| < 2i?. Hence 

< 



#1/6 |z — z„| 



giving 



5i < ci? 1/6 log(JJ) 



(32) 
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for some numerical constant c. In the sum S2 we use the same approach as before 
to obtain 

S* 2 < 16ei? 1/6 log(i?). (33) 
Combining (|3"Tj) . (|3^|) and (|3~3"1) we obtain 

< C 2 i?^ 1/6 + C*3Cxp(17cei? 1 / 6 )ei? 1 / 6 log(i?) 
in which C2 and C3 depend only on Q. We thus obtain the total estimate: 

\Kg(0,t) -K g (Q,t)\ < Cei? 1/6 log(i?)exp(17eei? 1/6 ) 

+ C(p-l)- 1 ^(l + \\q-q\\ p )mm(l, T L 

where C depends only on Q. Upon possible enlarging C equation ((9]) gives 

|fl(0,i)| < Cei? 1 / 6 log(i?)cxp(17cei? 1 / 6 ) 

logiT 
tR v 

The final result follows from LemmaEJwith d = CeR 1 / 6 log(i?) cxp(17 e eR 1 / 6 ) 
and Co = C(p — l) _1 / p (l + \\q — q\\ p ) min ^1, using the same method proof 



+ C(p-l)- 1 /*(l+\\q-q\\ p )mm( 



1. 



as for Theorem 

Theorem 6.1. Let Q\ and Q p be positive numbers and p £ (1,2]. Then there 
is a positive number C , depending only on Q\ and Q p , and a positive number Rq, 
depending on Q\, Q p , andp, so that the following statement is true for any R > Rq 
and any e € (0, 3/4). If q and q are two potentials in B(Q\) such that \\q — q\\ p < Q p 
and for which the respective zeros of the corresponding Jost functions are e-close in 
a disc of radius R, then 



sup 

x£[0,l] 



< 



CiQ^Q^^logR^-V^^R-^ 2 /^ 2 ^ 

+ ei?s log(i?) cxp(17 c eR^ ) | . 



Corollary 6.2. (Conditional stability). Let q and q be two potentials with support 
in [0,1]. Let J Q \q\ p < Q p and J Q \q\ p < Q p for some p > 1. Then for any 
8 > there exists a pair (e, R), depending only on 6, Q p , and p, such that if the 
corresponding Jost functions have zeros differing by at most e in a disc of radius R 
then 

sup / (q-q) <S. 
xe[o,i] Jx 

Theorem 16.11 and Remark [T^lso imply the next corollary. 

Corollary 6.3. (Uniqueness) Let q and q be two integrable potentials with support 
in [0, 1]. // all eigenvalues and resonances for one potential coincide with those of 
the other, then q equals q almost everywhere. 

Remark 3. Taking into account Remark[lJ one can prove conditional stability for 
q, q just from B(Q). But in this case the radius R can not be chosen uniformly. 
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Indeed, consider the sequence of potentials q n {x) := nX[o,i/n](x) where x is the 
characteristic function. Then for the corresponding Jost functions ip( z 'i n ) we have 



V z 2 — n \ . z . ( v 



2 



ib(z: n) = e %z l n [ cos I — ) — i sin 

K ' \ \ n J sfz^ 

Obviously, for any R there are no zeros of ifj(z; n) in the disc \z\ < R if n = n(R) 
is sufficiently large. Thus, one cannot choose R to depend only on S and Q. At 
the same time, in the same way as in the proof of the estimate (|27|) one can still 
get uniformity for the class of potentials from B(Q) if their behavior as x — » is 
specified. Namely, if we assume that q, q € B(Q) and 

(W\ + M) < V(e), 



o 



where 77(e) — > as e — * 0, then R can be chosen to depend only on 6, Q, and the 
function 77. 

Acknowledgment: The authors wish to thank Sergey Naboko for many discus- 
sions on the topic. 
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